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The spectral data on the diffuse Galactic 7-rays, at medium and high latitudes (| b |> 10°) and 
energies of 1-100 GeV, recently published by the Fermi Collaboration are used to produce a novel 
study on the 7-ray emissivity in the Galaxy. We focus on analyzing the properties of propagation of 
cosmic rays (CRs), using the publicly available DRAGON code. We critically address some of the 
models for the interstellar HI and H2 gas distributions commonly used in the literature, as well as 
test a variety of propagation models. Each model assumes a distinct global profile for the diffusion 
and the re-acceleration of CRs. Fitting propagation parameters to well measured local CRs such 
as, the B/C ratio, p, p and e ± fluxes, we evaluate the 7-ray spectra at medium and high latitudes 
in order to place further constraints on these propagation models. 



I. INTRODUCTION 

Interactions of CRs with the interstellar medium 
(ISM) are a copious source of 7-rays through 7r°, in- 
verse Compton and Bremsstrahlung emissions. The 
study of diffuse 7-rays at intermediate and high lati- 
tudes (|6| > 10°) d is a promising tool to probe and 
constrain the propagation of CRs in the Galaxy as well 
as ISM properties, with the CR propagation described 
by: 
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where tp(f,p,t) is the CR density, q(r,p,t) is the 
source term including components of primary origin, 
as well as CRs from spallation and decay processes 
from heavier elements. D xx (r) is the spatial diffusion, 
D pp (r) describes the diffusion in momentum space due 
to re-acceleration, p is the momentum loss rate due to 
interactions with ISM, the Galactic magnetic field or 
the interstellar radiation field (ISRF), V is the convec- 
tion velocity due to Galactic winds, while Tf rag and 
Tdecay are the timescales for, respectively, fragmenta- 
tion loss and radioactive decay. 

For our simulations we use the DRAGON code 0] 
that numerically solves Eq. ([1} in the steady state ap- 
proximation, in a 3D grid; 2 spatial dimensions for 
Galactocentric radial distance and height from the 
Galactic plane, and 1 for the momentum p. 



II. ASSUMPTIONS 



A. Primary Sources 



CR primary sources up to energies of ~ lOOTeV, 
are supernova remnants (SNRs). For each nucleus i 
the source term describing the injection of CRs in the 
ISM is given as a function of rigidity R by 



q l (r,z,E) = q ,if s (r,z)( 



Ro ' 



(2) 



where go,i is the normalization for each nucleus, 
f s (r\z) traces the distribution of SNRs as modeled 
in [3j on the basis of pulsar and progenitor star sur- 
veys [H . Electrons and positrons accelerated between 
a pulsar and the termination shock of the wind nebula 
may contribute to the high energy e spectrum, and 
then to the 7-ray flux, via inverse Compton scatter- 
ing. Following the parametrization of [5|, the source 
term due to a distribution of pulsars can be described 
by a power-law with an exponential cut-off given by 



Q p (r, z, t, E) = J E- n e- E / M f p (r, z) , 



(3) 



where Jo depends on the averaged birth rate of pulsars 
distribution and the average portion of pulsar initial 
rotational energy injected in the ISM as e ± , f p (r, z) 
describes the spatial distribution of young and middle 
age pulsars modeled by [6j, while n and M are, re- 
spectively, the injection index and statistical cut-off, 
for the pulsar distribution. 



B. Magnetic Field and Diffusion 

The large scale Galactic magnetic field is generally 
assumed to be a bi-symmetrical spiral with a small 
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pitch angle 0] • Here we assume that the regular mag- 
netic field is purely azimuthal, Bo = Bq4>, and has the 
form 



KfG), (4) 



based on the analysis of WMAP synchrotron intensity 
and polarization data in [8[. 

Diffusion is a result of CR scattering on randomly 
moving magneto-hydro-dynamical (MHD) waves. A 
correlation to the Galactic magnetic field and its dis- 
continuities is expected, with a larger diffusion coeffi- 
cient where the magnetic field is weaker. We assume, 
on the basis of 0), an isotropic diffusion coefficient of 
the form 



D(r,z,R) 



A)/3"(-^) 5 exp(- 

Kq 



exp( — ), (5) 

Zd 



where Rq — 3 GV is the reference rigidity, 8 is the 
diffusion spectral index depending on mechanism that 
builds up the turbulence in the ISM, and Zd are, 
respectively, the radial and vertical scales defining 
the diffusion profiles in the Galaxy. The dependence 
of diffusion on the particle velocity, f3 — v p /c, is 
naturally expected to be linear (77 = 1), however the 
analysis by [9( shows an increase in diffusion at low 
energies. To represent such a behavior, the parameter 
•q has been introduced by [To| . 



C. Interstellar Gas: 

The interstellar gas (ISG) is composed of hydro- 
gen, helium and small contributions from heavier ele- 
ments, with hydrogen observed in atomic (HI), molec- 
ular (H2) and ionized (HII) states. For the distribu- 
tion of HI gas we use as a reference the model devel- 
oped by 

Molecular hydrogen can exist only in dark cool 
clouds where it is protected against the ionizing stellar 
ultraviolet radiation. It can be traced with the A = 2.6 
mm (J = 1 — > 0) emission line of CO, since collisions 
between the CO and H2 molecules in the clouds are 
responsible for the excitation of CO. The CO to H2 
conversion factor, Xco which relates the H2 column 
density, Nh2 , to the velocity-integrated intensity of 
the CO line, has considerable uncertainties. For the 
H2 distribution we use for our reference model the 
map provided by [f2j , assuming the conversion factor 
to vary exponentially with Galactocentric radius, 



X C Q[H 2 cm 



^K^knT 



1 .4exp( 



I? 



(6) 



Tl(kpc)' 

however [l3| is also an older widely used model in the 
literature. The radial and vertical profiles of H2 distri- 
bution models are shown in Fig. [T] Ionized hydrogen 
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FIG. 1: Large scale molecular hydrogen distribution in the 
Galaxy vs r for z — O(left); vs z for r — r© (right), dashed, 
solid and dashed-dotted lines for our reference model [l^] 
with, respectively, high, mean and low values of H2 mid 
plane density, dotted lines for the model provided by 



occurs in the vicinity of young O and B stars, with 
the ultraviolet radiation from these stars ionizing the 
ISM. HII regions have a similar distribution to the 
molecular hydrogen, but mass-wise their contribution 
is negligible. 



D. Methodology 

The propagation parameters are determined upon 
fitting CR spectra. We consider a range of values for 
5, Zd and in Eq. For each set of these values 
we derive (Dq^AtIj) by minimizing the \ 2 for Boron 
over Carbon B/C data (Fig. [2] (upper left)). 

The injection spectrum of protons described by 
three spectral indices, "ff (see Table [J) is fitted to the 
PAMELA [3 and CREAM data [f| (Fig. [2] (upper 
right)). The predicted antiproton spectrum is consis- 
tent with local flux (Fig. [2] (lower left)). Helium spec- 
trum is also checked for consistency with the most 
recent data from PAMELA fiH ]. 

The e ± fluxes below E ~ 30 GeV is dominated by 
primary electrons accelerated by supernovae and by 
secondary electrons (and positrons) produced in in- 
elastic collisions of CR nuclei with the ISM. The spec- 
trum of the secondary e ± is related to the CR nuclei 
spectrum, while the primary electron spectrum can 
be described by a single power-law j e above 5 GeV; 
which wc fit from the low energy e + + e~ spectrum 
measured by Fermi (for a more detailed discussion see 
)■ Pulsars within ~ 3kpc may also contribute to 
the e + + e~ spectrum up to O(0. 1) at E 50 GeV 
and up to O(l) at E 500 GeV [f|. Assuming pulsars 
contribute maximally, we find the properties of pulsars 
Jo, n , and M in Eq. ([3]) upon fitting the Fermi data 
as shown in Fig. [5] (lower right). The positron frac- 
tion and electron spectra measured by PAMELA are 
checked for consistency as well. 
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FIG. 2: Upper left: B/C spectrum to fit the diffusion 
parameters in Eq. ([5]), Upper right: the proton flux to fit 
7? , Lower left: the predicted antiproton spectrum, Lower 
right: the e + + e~ flux to fit 7 e and pulsar parameters in 
Eq. ©. 




FIG. 3: Reference model, predictions for the 7-ray spec- 
trum for 0° < I < 360°, Upper left: 10° < \b\ < 20°, upper 
right: 20° < \b\ < 60°, lower: 60° < |6| < 90°. 



III. RESULTS 

Our reference model with S = 0.5, z c i — 4 kpc and 
rd = 20 kpc called "KRA4-20" provides a good com- 
bined fit of local CRs (Fig. [2]) and diffuse 7-ray spectra 
at intermediate and high latitudes (Fig. [3]) . 

In Fig. U we show the 7-ray spectra for differ- 
ent spectral indices; 8 = 0.5 ("KRA4-20"), 8 = 0.4 
("RUN4-20") and 5 = 0.33 ("KOL4-20"). The tt° 
component depends on the propagated proton spec- 
trum. Lower values of the 8 makes the propagated 
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FIG. 4: Predictions for different values of the diffusion 
index 8. dotted lines: S = 0.5, dashed lines: S — 0.4, 
dashed-dotted lines: 8 = 0.33. For all Zd = 4 kpc and 
rd — 20 kpc. 



FIG. 5: Predictions varying the diffusion radial scale rd- 
dotted lines: rd = 5 kpc, dashed lines: rd = 10 kpc, 
dashed-dotted lines: rd = 20 kpc. For all S — 0.5 and 
z d = 4 kpc. 



spectrum harder, resulting in the need for a softer in- 
jection index 7? in order to fit the data as shown in 
Table |U The differences in 7? + 6 causes small differ- 
ences in the 7r° fluxes. Electrons propagation, unlike 
protons, mainly depends on energy loss time-scale and 
thus is not affected much by varying the diffusion in- 
dex. Since we normalize our diffusion coefficient at 3 
GV (see Table UJ, for larger S the higher energy e dif- 
fuse faster, reaching the higher latitudes faster giving 
a slightly harder inverse Compton spectrum. Differ- 
ences in low energy Bremsstrahlung emissions come 
from different Alfven velocities (see Table HJ, with the 
greater re-acceleration depleting the low energy spec- 
trum. The decreasing of r d from 20 kpc ("KRA4-20") 
to 5 kpc ("KRA4-5") makes the diffusion coefficient 
smaller towards the Galactic center, which forces CRs 
produced closer to the Galactic center to spend greater 
time there. Since we refit the Do for each propaga- 
tion model (see Table HJ the net effect in the fluxes 
is negligible (Fig. [5]). In Fig. [6] we show the effect 
of varying the diffusion scale height z<i from 1 kpc 
("KRA1-20") to 10 kpc ("KRA10-20"). Since tt° and 
Bremsstrahlung emissions are morphologically corre- 
lated to the gas distribution which is concentrated 
close to the galactic disk, they do not change much 
by changing the size of diffusion zone. The inverse 
Compton spectrum is mainly affected by the actual 
distribution of electrons being confined within thinner 
(thicker) diffusion zones resulting in lower (higher) to- 
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FIG. 6: Predictions varying the diffusion scale height zj,. 
dotted lines: Zd — 1 kpc, dashed lines: Zd = 4 kpc, dashed- 
dotted lines: Zd = 10 kpc. For all 8 = 0.5 and = 20 kpc. 




FIG. 7: Predictions for different molecular hydrogen dis- 
tributions shown in Fig [I] 



tal inverse Compton flux, since even infrared and op- 
tical target photons have a much thicker distribution 
profile than the ISM gas. 

Among different components of the ISG, molecu- 
lar hydrogen distribution has large uncertainties (see 
Fig.[l}. Our reference for H2 distribution [12[ is called 
" NS mean" . To account for uncertainties in H2 mid- 
plane density we also study its high and low values 
called, respectively, "NS high" and "NS low". We 
study the model developed by [l3[ as well ("Bronf). 
The larger number of target nuclei in the gas model 
[IH , results in the need of a faster escape of CRs from 



Name 
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~ci 


TJ 


Do 


Va 


V 


ill ill ll 


R p i 


KRA4-20 


0.5 


4 


20 


2.49 


19.5 


-0.36 


2.06/2.35/2.18 


14.9 


KRA4-5 


0.5 


4 


5 


2.76 


16.9 


0.0 


2.07/2.35/2.18 


27 


KRA4-10 


0.5 


4 


10 


2.58 


19.1 


-0.25 


2.05/2.35/2.18 


17.5 


KRA1-20 


0.5 


1 


20 


0.55 


16.3 


-0.52 


2.07/2.34/2.18 


16.5 


KRA10-20 


0.5 


10 


20 


4.29 


19.1 


-0.37 


2.05/2.35/2.18 


15.2 


RUN4-20 


0.4 


4 


20 


3.21 


23.2 


0.32 


2.06/2.44/2.28 


14 


KOL4-20 


0.33 


4 


20 


3.85 


24.8 


0.77 


2.03/2.49/2.35 


10.7 



TABLE I: The parameters for the various models of prop- 
agation. See text, Zd and Td in kpc, Do in X 10 28 cm 2 s — 1 
(at 3 GV) and va in km/s; 7^ is the protons injection in- 
dex below the R\ (in GV), 7^ the injection index between 
R\ and R? 2 — 300 GV, and 73 above R%. For primary elec- 
trons we assumed one break at 5 GV above(below) which, 
the injection index is 2.62(1.6). The pulsar parameters are 
n = 1.4, M = 1.2 TeV. 

the Galaxy, in order not to overproduce secondaries. 
As a result, to keep the same B/C flux ratio, the dif- 
fusion coefficient normalization needs to be increased. 
Therefore e ± would propagate to larger distances from 
the Galactic disk, resulting in an increase in the ob- 
served flux. Thus the combined analysis of CRs and 
7-rays can constrain the large scale distribution of the 
gas in the Galaxy (see [r| for more details on ISM 
constraints). 

In conclusion, we have studied a rather extreme 
set of CR galactic diffusion profiles, ranging from an 
equivalently constant galactic diffusion coefficient, to 
a diffusion profile with large gradients both in r and 
z. Wc find that thicker diffusion zones are preferred 
by the combined fit of CRs and 7-rays, while the 
r-dependence of the diffusion coefficient can not get 
strongly constrained from the current data. More im- 
portantly, we find that even with the existing data, 
such a combined analysis can discriminate and even 
constrain profiles of the ISM gas. 
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